DFTT 26/96 
Cavendish-HEP-96/07 

May 1996 



°^ . 

2 ; Higgs signals and hard photons via VFVF-fusion 

in the Standard Model 
at the Next Linear Collider 



X 



Stefano Morettifl 



> 

00 

Dipartimento di Fisica Teorica, Universita di Torino, 

O . and I.N.F.N., Sezione di Torino, 

\0 ■ 

Q\ • Via Pietro Giuria 1, 10125 Torino, Italy. 



Cavendish Laboratory, University of Cambridge, 

P 



p Madingley Road, Cambridge, CBS OHE, United Kingdom. 



PACS numbers: 14.80.Bn, 14.70.Bh, 14.70.Fm, 29.17. +w. 



Abstract 

Within the framework of the Standard Model, integrated and differential distributions 
are given for Higgs production via the VbVF-fusion mechanism and decay via the channels 
H —* bb and H — > WW — * jjjj, with and without photon radiation, at Next Linear Collider 
energies. Calculations are carried out at tree-level and rates of the leading processes 
e + e~ — ► v e v e H — ► v e v e bb and e + e~ — > P e f e -ff — > v e v e WW — ► z^ejjjj are compared 
to those of the next-to-leading reactions e + e~ — > v e v e H{pi) — * VevJAr) and e + e~ — > 
v e v e H(^f) — > u e v e WW{^) — ► ^e^ejjjjT) i n the case of hard and detectable photons. Finally, 
a brief discussion concerning the case of H — > ZZ — ► jjjj (7) decays is also given. 



1 E-mails: Moretti@to.infn.it; Moretti@hep.phy.cam.ac.uk. 



1. Introduction 



The motivations for building an e + e~ linear collider operating in the energy range \/s = 300 — 
1000 GeV (NLC, Next Linear Collider) are indeed quite convincing |], @, [3|, |], [5| . For example, 
one of the major goals of such a machine will be to perform 'high precision' Higgs physics. That 
is, to measure the parameters (mass Mh, width Th, spin Sh, {J vc }h quantum numbers, etc.) 
and the relevant phenomeno logical rates (cross sections, branching ratios, etc.) of the Higgs 
boson H, to an accuracy that can not be achieved at hadron machines. The importance of this 
project is clear if one considers that the Higgs mechanism of spontaneous symmetry breaking 
is a sort of 'pedestal' that holds up the entire Standard Model (SAi). 

If it exists, the H particle will be probably discovered at the Large Hadron Collider (LHC) 
0, which is expected to start running early next century, well in advance of the future NLC. The 
CERN pp collider will however have great difficulties in measuring the fundamental quantities 
related to this particle, because of the huge background proceeding via strong interactions 
which will be present at the hadron machine. This is particularly true for a Higgs boson in the 
intermediate mass range (IMR), Mh ~ 2Mjy, for which the detectability of the only two viable 
decay channels, H — > 77 and H — > bb, strongly depends on the detector performances. For a 
Higgs in the heavy mass range (HMR), Mh ~ 2Mw, things should be easier, as the four-lepton 
decay mode H — > ZZ — > A£ is relatively clean and straightforward to detect. 

At the NLC, the main SAi Higgs production mechanisms are via the bremsstrahlung process 
e + e~ — > ZH and via the WW- and ZZ-fusion reactions e + e _ — > v e v e WW (e + e~ Z Z)^> 
£/ e i/ e (e 4 ~e~)H ||. For a first stage NLC (with y/s ~ 500 GeV) the rates of the bremsstrahlung 
mechanism are larger than those of the fusion channels, if Mh ~ 2Myy. At larger centre-of- 
mass (CM) energies (y/s ~ 500 GeV) H^H^-fusion starts dominating over the whole of the Mh 
spectrum. Furthermore, the rates of ZZ-fusion are generally one order of magnitude smaller 
than those of the IViy-channel. Concerning possible Higgs signatures, all the principal decay 
modes of this particle can be detected and studied at the NLC, for all values of Mh, provided 
that enough statistics can be accumulated . In particular, the main sources of Higgs events 
will be the channels H -> bb (IMR) and H -> WW -> 4jet (HMR) 0. 

There are a few important aspects, in the way a NLC running around the TeV energy scale 

will operate, that should be carefully considered in order to exploit in full its potential and 

that are absent in lower energy e + e~ machines. These are related to the influence on the cross 

sections of electromagnetic (EM) interactions which can take place before the actual beam 

collision, such as bremsstrahlung and beamsstrahlung effects [00. Whereas effects due to 

synchrotron radiation emitted by one of the colliding bunches in the field of the other one (i.e., 

2 Effects due to the energy spread of the beams before annihilation (intrinsic to any collider) should also be 
considered. 
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beamsstrahlung) necessarily need, in order to be quantified, the knowledge of the technical 
details of the collider design and can be realistically estimated only through Monte Carlo 
simulations, those due to the emission and exchange of photons from and between the actual 
pair of electron and positron which collide (i.e., bremsstrahlung or Initial State Radiation, 
ISR), can be treated in a fairly general way. In many cases, one can compute the exact EM 
corrections to the e + e~ annihilation subprocess and express these via the so-called 'electron 
structure functions'. Such corrections embody both real and virtual photon radiation and they 
are known to date up to the order 0(a^ m ) |12]. It has also been shown that, for 'narrow beam' 



designs, beamsstrahlung affects the cross section much less than the ISR ||11|| , such that in 
phenomenological analyses one can consistently confine oneself to dealing with bremsstrahlung 
radiation only. In general, the principle effect of the ISR is to lower the effective CM energy 
available in the main process, thus ultimately reducing (enhancing) total cross sections which 
increase(decrease) at larger CM energies. Furthermore, ISR also leads to a smearing of the 



differential distributions |13|. However, the 'electron structure functions' approach is not always 
applicable, such as in the context of H^VT-fusion processes in e + e~ annihilations. 

It is the purpose of this paper to study the properties and the effect on the integrated rates 
as well as on the differential distributions of interest to Higgs searches of hard photon emission, 
which can take place in the SM. Higgs production and decay processes 

e + e" -> v e v e H(j) -> v e v e W), (1) 

e+e~ -> v e v e E{^) -> p e v e WW(j) -> ^ejjjjT, (2) 

(that is, via the lyiy-fusion mechanism), at NLC energies and for several values of the Higgs 
mass, and where the photon can be emitted from the initial, virtual and final states. In our 
opinion, the importance of this study is motivated by the fact that, contrary to lower energy 
e + e~ colliders (such as LEP1, SLC and in part also LEP2), which 'sit' on gauge boson resonances 
and so the width of the unstable particles imposes a natural cut-off on events with hard photons 
produced by the initial state, at the NLC such a suppression does not act any longer. In 
addition, as the beam energy is much larger, the probability that the incoming electrons and 
positrons can radiate hard photons increases greatly. From this, it follows that in practise a 
sample of pure e + e~ — ► v e v e bb and e + e~ — > Pe^ejjjj events, without 7-radiation, does not exist 
and one inevitably has to deal with EM emission. Moreover, in the experimental samples of 
data, events in which the photon is emitted during the production process e + e~ — > v e v e H^ 
are not distinguishable from those in which the radiation comes from the Higgs decay stages 
H — > bby and H — > WW^) — *> jjjj7, such that in phenomenological simulations one is forced 
to consider both the cases at the same time. In particular, we stress that it is principally the 
second kind of radiation which could spoil the shape of the Higgs resonances. 
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The plan of the paper is as follows. In Sec. 2 we devote some space to illustrate the 
computational techniques that we have used as well as the numerical values adopted. Sec. 3 
presents our results, and in Sec. 4 we give a summary and report the main conclusions to be 
drawn from this study. 



2. Calculation 

The Feynman diagrams corresponding to processes (|l])-(@) are given in Figs. 1 and 2, respec- 
tively. To compute them we have adopted numerical techniques which use the helicity amplitude 
formulae of Ref. [jH and the subroutines contained in the package HELAS [T5|] . For process 



([!]) we have also resorted to the program MadGraph [16| for the generation of the FORTRAN 



code of the Matrix Element (ME). The routine VEGAS [17] has been used as multi-dimensional 
integrator. To achieve high accuracy in the evaluation of the cross sections in presence of Higgs 
peaks in different regions of the phase space, we have used the method of splitting the Feyn- 
man amplitudes squared into non-gauge invariant terms, each of which with a different peak 
structure. These have been then integrated separately by using an appropriate choice of the 
phase space variables, according to the resonant behaviour of the graphs involved. A sum over 
the various terms of the MEs gives in the end gauge-invariant results. Such a procedure has 
been already described elsewhere (for example, in Ref. |18|), such that we do not enter here 
into details. 

The computation of the graphs in Figs. 1-2 is rather straightforward, though the numerical 
integrations over the phase space can in some instances be rather complicated, as they can 
involve up to 16 dimensions^]. Thus, in order to perform our calculations in a reasonable amount 
of CPU time we have always constrained the iy-bosons to be on-shell in process (0)0. In this 
way, we are able to reduce by two the number of the variables we integrated over. As we will be 
eventually interested in studying the invariant mass distributions of the decay products of the 
Higgs boson, we expect such a procedure to have a little impact on our final results. However, a 
direct consequence of this simplification is that we are not able to study Higgs decays into four 
jets below and near the 2M\y threshold. Nonetheless, since the phenomenology of the off-shell 
Higgs decay H — > W*W* is well known, predictions in the range M# ~ 2Mw can be easily 
extrapolated from our rates. For consistency, we have constrained the ly-bosons to be on-shell 



3 Corresponding to a final state of 7 particles, in case of process (| 
4 Both in the case of process (0) and in the lowest order reaction 



e + e — * v e v e H — > v e y e bb 1 (3) 
no simplification in any respect has been adopted. 
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also in the case of the leading process 

e + e~ -> v e v e E -> v e v e WW -> z/ e i/ e jjjj. (4) 

Furthermore, the interferences between the various peaks in processes (|l|)-(@) are in general 
very small compared to the squared terms of the resonances and do not bring any distinctive 
structure into the differential distributions. In fact, such contributions always mix up graphs 
with different resonant structures^ such that the phase space regions in which one or more 
amplitudes are large are different for different graphs. Therefore, as a further simplification, 
such interferences have been systematically neglected in the presentation of our results. 

The following numerical values of the parameters have been adopted: Mz = 91.175 GeV, 
Tz = 2.5 GeV, Mw = 80.23 GeV, and for the Weinberg angle we have used its leptonic effective 
value sml s (8w) = 0.2320. For the fermions: m e = m Ue =0, m;, = 4.25 GeV, whereas all light 
quarks u, d, s and c have been considered massless. Jets have been identified with the partons 
from which they originate. The EM coupling constant a em has been set equal to 1/128. For 
the Higgs width Th we have adopted the tree-level expression corrected for the running of 
the quark masses in the vertices Hqq (these have been evaluated at the scale \i = M H fl9fl ). 
Therefore, in order to be consistent, we have used a running 6-mass in the H — > bb vertex of 
the production processes ([!]) and (|3|). As representative values of the CM energy of the NLC 
we have adopted 300, 500 and 1000 GeV, and the Higgs mass has been chosen in the range 
60 GeV ~ Mh ~ 450 GeV. 

3. Results 

Photon radiation can be emitted in processes ([!])— (^) both during the production mecha- 
nism e + e~ — > v e v e H^i and during the decays H — > bby and H — > ^ e ^ej j j j 7 • The emission 
is gauge-invariant in both these two stages and we will refer to it as 'production radiation' and 
'decay radiation', respectively. The ISR as previously defined would correspond here to the 
photons emitted by the incoming electron and positron lines. The diagrams describing such 
bremsstrahlung (i.e., graphs 1-2 in both Fig. 1-2) do not satisfy gauge- invariance if taken on 
their own. To recover the latter it is needed to add the graphs in which the photon is emitted 
by the W^-lines (i.e., graphs 5-6 in both Fig. 1-2). Therefore, a separation of the ISR diagrams 
from the rest does not make sense here. The issue of defining the ISR in a gauge-invariant 
way in presence of charged current (CC) interactions is indeed a well known problem, for ex- 
ample in case of e + e~ — > W + W~ 4 fermions production at LEP2. In the sense that, in 

5 We have integrated them by using a flat phase space, which does not map any of the possible peaks of the 
interfering graphs. 
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presence of electric charge flow between the initial and the subsequent stages of e + e~-initiated 
processes, the definition of the ISR is not unique. Such a difficulty has been probably overcome 
in the case of neutrino exchange, by means of the so-called current splitting technique (for 



details, see Ref. p0|). In that case, complete ISR separates into a universal, factorising, process 
independent-contribution (that is, the electron structure functions of Ref. [T2|]) and a non- 
universal, non-factorising, process-dependent part (that is, the complex analytical expressions 
as given in Ref. |21|). However, this approach is not realistically possible in the present context. 



In fact, the formulae of Refs. |I2, are valid only in case of annihilation- and conversion- type 
Feynman diagrams, that is, when the incoming e + - and e~-lines are connected to each other 
via a s- or t, w-channel, respectively. In contrast, we are concerned here with electron/positron 
lines that are disconnected from each other and that end up as the (anti)neutrino ones of the 
final states, by means of multiple space-like CC-interactions. 

The most correct approach would certainly be to compute the complete (D(a em ) (and be- 
yond) corrections, including the loop diagrams. That is however beyond the intentions of this 
study. What we will do here is to show results for the leading order (LO) processes (0)-(f|) and 
the next-to-leading order (NLO) ones (|l])-(@) separately, the latter with real photons having 
Pr > 1 GeV (hard and detectable EM radiation), in terms of both integrated and differential 
rates. We also impose an additional cut on the photon, by asking that the latter be isolated 
from the jets arising in the final states of processes (H])-©, from the b- and the light quarks: for 
example, by imposing cos 6^ - n < 0.95 (corresponding to a cone with an angular size of ~ 18 
degrees). This is done because it would generally be impossible to tag photons too close to the 
original parton, as the latter gives rise to a jet with a finite angular size, such that if the photon 
fails within the corresponding cone it will not be distinguished from the other parts of the jets. 
In this case, its energy is counted as part of that of the hadronic system associated with the par- 
ton and the M b i and Mjjjj invariant masses are not experimentally measurable in events of the 
type (|l|)-(@). For configurations in which the photon is highly collinear (i.e., cos 6^ j 7 > 0.95), 
such events would rather be recognised as leading 2^4 and 2^6 topologies. 

In order to be consistent with the sketched approach, we will refrain from summing up lowest 
and next-to-lowest rates through the orders 0(a^. m ) and 0(al m ), as such a summation would 
need to include also the contributions due to virtual photons. In practise then, the aim of our 
study is, on the one hand, to estimate the size of the rates produced by events with hard photons 
and, on the other hand, to compare their kinematic properties to those of the non-radiative 
processes. This is done in order to establish whether the shape of the differential spectra of 
interest to Higgs searches can be significantly modified at higher order by the presence of hard 
photons, such that when proceeding to experimental analyses one might expect a smearing of the 
relevant resonant distributions. In particular, we remind the reader that the contributions that 
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we have not computed (the loop diagrams) or removed by the cuts (the infrared photon regions) 
would show the same kinematics as the lowest order processes, such that they would modify 
the overall normalisation of the differential distributions but not the shape. In particular, 
concerning the 'decay radiation', it is well known that the Kinoshita-Lee-Nauenberg theorem 
Z2| | prescribes that all the logarithmic contributions of the form (a em /7r) n ln n (s/m^ ■), which 



appear to the order 0{a em ) n according to the Sudakov theorem [gjj and which are due to 
'collinear' emission of photons from the final state, must not appear in the expression of the 
inclusive cross section. They are in fact canceled by the negative contributions due to virtual 
final state photons. Therefore, the hard photon effects that we will discuss in the following 
could well be larger in the complete 0(ot em ) result. In the very end, in order to compare theory 
and experiment, the latter are certainly needed. However, our preliminary results will enable 
us to assess whether, at higher order, complications should be expected in individuating the 
position of the peaks and/or in establishing their line-shape, or whether the effect of 0(a em ) 
corrections is mainly matter of overall normalisation of the leading order distributions. 

In our analysis, we will present spectra in M 6 g, M fe g 7 , Mjjjj, (that is, in the invariant 

masses of the decay products bb(j) and jjjj (t) °f the Higgs scalar), together with those in E 1 
(the energy of the photon) and in p™ lss (the missing transverse momentum of the undetected 
neutrino pair). The importance of the distributions in invariant mass is clear if one considers 
that the presence of 'decay radiation' could spoil in the end the form of the Breit-Wigner peaks 
in M b i and Mjjjj that one obtains from the LO contributions (|3])-(fD, via diagrams 3-4 in Fig. 1 
and 3-4 & 7-10 in Fig. 2. On the other hand, these diagrams contribute to produce Higgs 
resonances in the M b ^ and Mjjjj 7 spectra. In contrast, the EM emission that takes place via the 
'production radiation' should not significantly spoil the form of the peaksf] and should instead 
contribute to enhance the number of resonant events. The distributions in E 1 and p™ lss could 
in principle be useful in order to separate in events of the type contributions due to 

'production radiation' from those due to 'decay radiation', provided significant differences exist 
between the corresponding differential spectra. In particular, these can be exploited in order 
to enrich the experimental sample of events with photons not emitted during the Higgs decays. 
Also, their knowledge could eventually be useful in disentangling Higgs signal from background 
processes (both irreducible and reducible ones). 

Finally, the cut on the transverse momentum of the photon has actually also been applied 
to all other particles in the final states of reactions ([j])— @j, apart from the neutrinos. Such 
a constraint should generally meet the requirements due to the finite coverage in energy and 
polar angle of the NLC detectors. 



^Its effect being mainly confined to lower the invariant mass of the W^VF-annihilation into the Higgs boson. 
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3.1. The Higgs process e + e — > v e v e blq 

We first consider the bb Higgs decay channel. The cross sections for processes (p]) and (|3[) are 
given in Fig. 3, for the mentioned choice of CM energies. They are presented as a function 
of the Higgs mass Mjj. For illustrative purposes we have considered the mass range 60 GeV 
~ Mh ~ 350 GeV, however the largest rates occur for Mh ~ 140 GeV (intermediate mass range), 
where the H —> bb decay has the largest branching ratio (BR). For M H > 140 GeV the off-shell 
channel H — ► W*W* starts dominating the SM. Higgs decay phenomenology. This is reflected 
in the steep decrease of the cross sections, both at leading and next-to-leading order, slightly 
before the real 2Myy threshold. 

The main feature of Fig. 3 is certainly the relatively large value of the rates for the radiative 
process ([!]) compared to those of the non-radiative reaction (^). In fact, the former is at least 
10% of the latter over all the interesting range of Mh, at all values of y/s. It increases with 
the CM energy and it is largely independent of the Higgs mass. In particular, if one assumes 
an integrated luminosity of / Cdt = 10 fb _1 , more than 300 events with hard photons could be 
produced per year, at a NLC with \/s = 1000 GeV and for M H rs 60 GeV. Such a number gets 
smaller with decreasing CM energy and increasing Higgs mass. 

The kinematic properties of events of the type ([!]) and (|) are illustrated in Figs. 5a-c, 
for yfs = 300, 500 and 1000 GeV, respectively. The selection of Higgs masses we have chosen 
here is Mh = 60, 100 and 140 GeV. For these values, the Higgs width is rather small, Th ~ 8 
MeV, such that hard photon emission in the decay process H — ► 667 (diagrams 3-4 in Fig. 1) 
is heavily suppressed. In fact, the contribution due to the two mentioned diagrams to the total 
cross section in e + e~ — > v e v e H{^f) — > v e v e bb^f events is rather small (by more than one order 
of magnitude, for all y/s and M H combinations) compared to the contribution due to graphs 
1-2 & 5-6 in Fig. 1, that is when the photon is emitted by the electron/positron and gauge 
boson lines. This is clearly reflected in the NLO M b i invariant mass distribution (upper left in 
Figs. 5a-c). In fact, the smearing towards low masses of the Higgs peaks due to the photon 
emission in the radiative H — ► bbj decay is essentially irrelevant for realistic phenomenological 
analyses: at the level of O(10~ 2 ) or even lessQ. 

The suppression of the 'decay radiation' is also visible in the invariant mass of the bbj 
system (upper right plots in Figs. 5a-c). In the sense that the corresponding spectra show 
a step (rather than a peak) at M b i y fa Mh, with a long and quantitatively relevant tail for 
M bbl > Mh, due to the photons from the 'production radiation'. In practise, the latter largely 
'overwhelm' those originating in the 'decay radiation', which would produce a Breit-Wigner 
peak. 

The distribution in energy of the radiated photon is displayed in the bottom left plots of 
7 Note the logarithmic scale in all the figures. 
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the same figures. It is qualitatively the same regardless of the actual value of the Higgs mass. 
In particular, very hard photons are suppressed at smaller (see Fig. 5a) and enhanced at larger 
(see Fig. 5c) CM energies, while the dependence of the spectrum at 500 GeV (see Fig. 5b) is 
roughly exponential. An additional suppression comes with the increase of the Higgs mass, 
at fixed y/s, especially if y/s < 500 GeV, whereas for y/s = 1000 GeV such an effect is quite 
small. Furthermore, as one expects, the 'decay radiation' photons are always softer than the 
'production radiation' ones. 

The distribution in missing transverse momentum is rather similar for both processes (|I|) 
and (HD (bottom right of the mentioned figures) and for the two radiative components at NLO 
(for this reason the latter are not shown separately here). In general, the p™ lss spectrum is 
only slightly harder for the non-radiative process. The maximum of the distribution is more 
pronounced as the CM energy increases and does not depend on the values of Mh- 

3.2. The Higgs process e + e~ — > v e Ve))]]l 

Integrated cross sections at the CM energies y/s = 300, 500 and 1000 GeV for processes (0) and 

@ are shown in Fig. 4, as a function of the Higgs mass in the heavy mass range|. Whereas at 

lower values of y/s (upper plot) the production of a Higgs scalar with Mh ~ 180 GeV is heavily 

suppressed, at larger CM energies (central and lower plots) rates are high enough to produce 

Higgs bosons up to masses of 450 GeV or so. In particular, at y/s = 1000 GeV, cross sections 

for both processes (|2|) and (|4]) diminish by only a factor of 8 if Mh increases from 180 to 450 

GeV. In contrast, at y/s = 300 GeV, there is a steep decrease of the rates as Mh approaches y/s. 

Curiously, in this case, the shape of the curves is rather similar to those of processes (|l|) and (§) 

(Fig. 3 upper plot). However, there the effect was due to the Higgs branching ratio, whereas 

here it occurs because of a phase space suppression on the production mechanism. Again, the 

most remarkable feature of Fig. 4 is a large rate for the radiative reaction (0), compared to 

the lowest order one (|J). This is generally true for all relevant combinations of masses and 

energies. The ratio between the rates of the two processes is less than a factor of 10, and up to 

100 radiative events per year can be expected (lower plot for the minimum Higgs mass). 

Figs. 6a-c show the kinematic properties of the particles in the finale state of reactions (0J 

and (f|). Since in the heavy mass range the Higgs width is quite large (it varies from « 0.6 

GeV at 180 GeV to ~ 40 GeV at 460 GeV !), hard photon radiation in the Higgs decay is 

no longer suppressed (unlike the case of H — > btry radiative decays), especially at large Higgs 

masses. Moreover, in process (§), there are six graphs which contribute to produce a peak in the 

8 We consider values of Mh (greater than 180 GeV) far above the real WW threshold at 2M\y ~ 160 GeV, 
in order to avoid complications due to having adopted here a narrow width approximation for the W boson, 
when the real value of Tw is instead around 2 GeV. 
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invariant mass of the jjjj7 system (i.e., number 3-4 & 7-10 in Fig. 2), since in the decay process 
H — > WW — > (7)jjjj7 photon emission can take place off two boson and four fermion lines. 
These two aspects (i.e., large Higgs width and 'multiple' emission) act in such a way that in the 
end the amplitude squared corresponding to the contribution of the six mentioned diagrams 
('decay radiation') is comparable to that of the others graphs ('production radiation'), over all 
the Mn spectrum and for all values of y/s considered here. This can be clearly appreciated by 
looking at the spectra (Figs. 6a-c, upper right). All the distributions in invariant mass 

of the jjjj7 system show a clear peak around the selected values of Mh (upper right plots in 
Figs. 6a-c). When the Higgs boson is lighter, such that its width is relatively narrow (for 180 
GeV in the plots), many hard photons still come from the 'production radiation', giving the 
long tail especially visible in the continuous lines, similarly to process (|1|). However, contrary 
to the case of H — > bb'j resonant decays, here the peaks at « Mh remain visible. When 

instead the Higgs width is larger (for Mh ~ 220 GeV in the plots), a large part of the photons 
come from the 'decay radiation', such that they contribute to build up a rather clear resonant 
dependence and the tail at M^ > M H starts disappearing (dotted and dashed lines). It is 
worth noticing that a phase space effect also contributes to enhance the above suppression at 
large values of Mjjjj 7 , such that in the end, in some instances (for example, at y/s = 300 GeV 
and Mh ~ 220 GeV, Fig. 6a), the NLO Mjjjj 7 and Mjjjj distributions look pretty similar. 

When considering the case of the distributions in Mjjjj at NLO, effects due to phase space 
enhancement /suppression play a more determinant role. In the sense that, when the portion 
of phase space available is restricted (that is at y/s = 300 GeV, Fig. 6a upper left, the shaded 
curves) the main effect due to the hard photons is the appearance of a tail at Mjjjj < M# 
(see especially the cases M H = 220 and 260 GeV). This is a consequence of a superposition 
of a rather symmetrical resonance (due to the 'production radiation' diagrams, whose shape 
looks similar to the one at LO, appropriately rescaled) and of a spectrum (due to the 'decay 
radiation' diagrams) that is strongly shifted towards low invariant masses^. As y/s increases 
the low mass tail effect due to diagrams 3-4 & 7-10 in Fig. 2 is counterbalanced by the fact 
that the contributions at NLO to the Mjjjj spectra due to diagrams 1-2 and 5-6 in Fig 2 have 
no longer a symmetrical shape (compare to the same distributions at LO, especially in Fig. 6c), 
but this is significantly shifted towards high masses. In many cases, in fact, the phase space 
available for Mjjjj > Mh is much larger than that in the complementary region Mjjjj < Mh- 
Therefore, for y/s ~ 500 GeV (Figs. 6b-c), the overall effect is a rather symmetrical resonance 
at NLO too. The same effect was not quantitatively appreciable in the case of process ([!]), since 
there the small width of the Higgs boson in the intermediate mass range drastically reduced 

9 Note that the appearance of a low mass 'tail' also for the LO distribution in Mjjjj when Mjj = 260 GeV at 
y/s — 300 GeV (upper dashed line in Fig. 6a, upper left plot) is accidental, as in this case the spectrum is close 
to its upper kinematic boundaries. 
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the relevance of configurations in which M b i — Mh 3> r# . 

The energy spectra of the emitted photon are shown in the bottom left plots of Fig. 6a-c. 
Again, one can see that very hard photons are not allowed if the Higgs mass/width and the 
CM energy are small (continuous line in Fig. 6a) and that photons created by the 'production 
radiation' are in most cases harder than those generated by the 'decay radiation', although the 
difference is less evident than for process (|l|). Unfortunately, this is not true at yfs = 300 GeV, 
where the distribution of the full sample closely resembles in shape that produced by diagrams 
3-4 & 7-10 in Fig. 2. 

The behaviour of the spectra in p™ lss (bottom right frames of Fig. 6a-c) shows that, also 
in the case of processes (0) and (f|), the average value of the missing transverse momentum is 
approximately the same both at leading and at next-to-leading order. Moreover, if one studies 
separately the contributions to the NLO spectra due to the 'production radiation' and 'decay 
radiation' diagrams, one realises that once again they do not show any substantial difference, 
their shape being described by the lower curves in Fig. 6a-c, with appropriate normalisations. 

Finally, a last comment is in order if one considers that a heavy H boson can decay into 
four jets also via the ZZ channel. Since BR(W — > jj) ~ BR(Z — > jj) (when also 6-jets are 
considered in the Z decay), the contribution to the total rates of the 2 — > 6 leading order 
reaction e + e~ — > v e i> e H — > i> e i> e \WW + ZZ] — > z/ e z/ e jjjj due to diagrams involving the H — > ZZ 
decay is proportional to the corresponding branching ratio. In particular, in the heavy mass 
range investigated here, these are expected to be smaller by a factor of ~ 20, for Mh = 180 
GeV, and of ~ 2.5, for Mh ~ 250 GeV, compared to rates due to contributions proceeding 
through the H — > WW channel. At next-lo-leading order, that is for the 2^7 reaction 
e + e~ — > v e v e H{^{) — > u e i> e \WW + ZZ](j) — > t'e^ejjjjT, a further suppression occurs. In fact, 
when H — > ZZ — > jjjj7, diagrams 3-4 in Fig. 2 no longer appear in the matrix element. 
However, the contributions of the corresponding squared amplitudes to the total cross sections 
of the above process are smaller than those of graphs 1-2 & 5-6, as well as than those of graphs 
7-10 (Fig. 2). It varies from a few percent at small energies (i.e., y/s = 300 and 500 GeV) and 
small Higgs masses, up to ~ 15% for yfs = 1000 GeV and large values of M#[^]. Thus, in the 
end, the inclusion into the calculations of diagrams involving H — > ZZ j j j j (t) decays could 
in some instances reduces the relative importance of events with hard photons. Conversely, this 
also indicates that experimental procedures exploiting the H — > jjjj decay inclusively, or the 
H — > ZZ — > jjjj on its own, are more likely to be less sensitive to higher order electromagnetic 
effects. 

10 We also expect an additional suppression for diagrams 7-10 of Fig. 2 when Mh ~ 2Mz- 
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4. Summary and conclusions 



In this paper we have studied the two radiative Higgs processes 

e + e~ -> v e v e H{i) -> u e u e bb-f, 

e + e~ -> v e v e H("{) -> v e v e WW{i) -> v e Ve}}}h- 

We have considered the case of hard and detectable photons, by adopting the cuts j ' 7 > 1 
GeV, on the transverse momentum of all particles in the final states, and cos 0^ j7 < 0.95, 
for the separation between the photon and the hadronic systems. We have then compared 
integrated and differential rates obtained for the two above processes to those of the reactions 
{p h i j > 1 GeV) 

e + e~ — > v e v e H — > z/ e z/ e 66, 

e + e~ -> z> e z/ e # -> v e v e WW -> v e v e }]}]- 

These four processes represent the main tree-level production and decay mechanisms of a Stan- 
dard Model Higgs boson in the M^VT-fusion channel, at the centre-of-mass energies typical of 
an e + e~ Next Linear Collider, at next-to-leading and leading order in the electromagnetic cou- 
pling constant, respectively. The final states v e v e bb(pj) have been considered in order to give 
account of the most likely decay mechanism of the Higgs boson in the intermediate mass range 
Mh ~ 2Mw (into a bb pair), whereas the ^ e ^e]]]]{l) ones refer to the favourite Higgs signature 
in the heavy mass range M ~ 2Mw (into four jets via WW). The values 300, 500 and 1000 
GeV have been considered for y/s. 

Our study has been motivated by noticing that the presence of hard electromagnetic emis- 
sion via bremsstrahlung photons from the initial state is somewhat unavoidable at these ma- 
chines. In fact, whereas at current lower energy leptonic colliders (such as LEP1 and SLC, 
and partially LEP2 as well) the width of the unstable particles produced in the e + e~ direct 
annihilation subprocess (Z and WW, respectively) naturally imposes a cut-off on hard photon 
radiation from the incoming electron and positron lines, this is no longer the case at the Next 
Linear Collider. Therefore one inevitably has to deal with such electromagnetic effects when 
attempting phenomenological analyses. Furthermore, in experimental data samples, radiative 
events in which the photon comes from the Higgs production mechanism are not separable from 
those in which it comes from the Higgs decay channels. Whereas in the former the shape of the 
Higgs resonances is generally not spoiled by hard photons, in the second this can in principle 
be heavily distorted. 

Since in our computation the contributions due to the loop diagrams were missing, we 
have not performed a summation through the orders 0(a^ m ) and 0{a 7 em ). Certainly, definite 
estimates of the 0(a em ) effects should be given only after a complete calculation including 



11 



virtual photons contributions. However, the following interesting indications can be extracted 
from our partial results. 

• Rates at next-to-leading order involving hard photons can be up to 30% of those at leading 
order. In particular, for all the relevant Higgs masses and CM energies considered here, 
these are never smaller than 10%. This corresponds to a rather large effect, in some 
cases comparable to or even larger than those due to beam related phenomena (such 
as beamsstrahlung, Linac energy spread), which are expected to take place at the high 
energy e + e~ colliders of the next generation. 

• The Higgs width acts as a regulator of the size of the EM emission in the Higgs decay 
processes. The smaller it is, the more suppressed the contribution to the total rates of 
diagrams involving photon radiation after Higgs production. Therefore, the smearing of 
the resonant distributions can be quantitatively significant only in case of Higgs masses 
in the heavy range (i.e., M H £ 2M W ). This should remain true also for poor detector 
resolutions. 

• If one considers heavy Higgs bosons, then two different scenarios appear. 

1. For a NLC with y/s ~ 500 GeV, the higher order distributions relevant to Higgs searches 
and studies (i.e., the invariant masses of the Higgs decays products bb and jjjj) generally 
have a shape similar to that of the 2^4 and 2^6 leading processes. This is because 
the smearing of the Higgs peaks towards low masses due to the emission of hard photons 
in the Higgs decays is counterbalanced by a tail towards large masses due to the large 
portion of phase space available to photons emitted in the Higgs production mechanism. 
In this case, the inclusion of the complete 0(a ern ) corrections is merely a matter of a 
different normalisation of the lowest order distributions. 

2. In contrast, for a \/s = 300 GeV NLC, where such phase space effects are negligible, the 
resonant differential spectra at NLO do show a quantitatively significant low mass tail. 
The effect clearly increases with the Higgs mass/width. In this case, as rates of events 
with hard photons can be large, a complete 0(a em ) calculation is desirable in order to 
assess the correct shape of the distributions at the complete NLO. 

• Since events with hard and detectable photons are defined on their own once appropriate 
cuts are applied to remove the infrared divergences (and these can well coincide with those 
dictated by the tagging procedures of the experimental analyses), it is feasible to study 
such events separately. An interesting features of these is that the invariant masses of the 
systems blry (for intermediate mass Higgses) and especially jjjj7 (for heavy mass Higgses) 
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show a clear step/peak around the actual value of M#. In many cases, a detectable 
number of events in the resonant region can be produced after a few years of running at 
the nominal collider luminosity. 

• Distributions which could naturally allow a separation between radiative events with 
photons emitted in the production (on the one hand) and decay (on the other hand) stages 
(such as, for example, the energy of the photon and the transverse missing momentum 
of the neutrino pair) are not generally helpful. The p™ lss spectrum presents the same 
behaviour for all energies and Higgs masses, both for the leading process and the two 
next-to-leading radiative contributions. The _E 7 spectrum could be exploited in most 
cases, but unfortunately not when y/s = 300 GeV and M# ~ 2M^, since in this case the 
spectra of the 'production radiation' and of the 'decay radiation' are extremely similar. 

In general, we stress that, in order to counterbalance the positive infrared divergences of soft 
and collinear (to both electron/positron and partons in the initial and final states, respectively) 
real EM emission, the corrections due to virtual photons in loop diagrams should generally be 
negative and could well reduce the relative contributions of events with the kinematics of the 
lowest-order processes. Therefore, the hard photon effects predicted here could well be enlarged 
by the complete result. 

In this paper, we have also pointed out that in four-jet decays of heavy Higgs bosons 
proceeding through the ZZ channel (i.e., H — ► ZZ — > jjjj7) hard photon effects are smaller 
than in the H — > WW^) — > jjjj7 case, but only over appropriate regions of the M H range 
(that is, well above the 2Mz Higgs decay threshold). 

Finally, some technical details. Our computations have been performed by resorting to 

a numerical evaluation of the exact matrix elements of all the above processes, by means 

of helicity amplitude techniques, integrating them over the appropriate phase spaces, using 

a standard Monte Carlo routine. A dedicated treatment of the various resonant diagrams 

makes the FORTRAN code produced rather accurate and fast. In this respect, we emphasise 

that we have used no approximations here, apart from that of constraining the two W bosons 

entering in the Higgs decay channel H — > WW^) — > j j j j (t) to be on-shell. This has been 

done in order to reduce the total amount of CPU usage, especially needed when performing 

a N = 3n — 5 (assuming no transverse polarization of the beams) dimensional integration 

for an n-particle final state. However, the knowledge of the off-shell decay phenomenology of 

the Higgs boson into two W's makes straightforward the extrapolation of our results to the 

relevant Mh rangeQ. Furthermore, some of the interferences between Feynman diagrams with 

11 Alternatively, the constraints on the invariant masses flowing through the W-resonances can be removed, 
and the integration performed over the original N dimensions, at the cost of additional CPU time. 
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different resonant structures have been neglected in the presentation of our results, as these are 
numerically negligible for the values of y^s and M H adopted here. 
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Figure Captions 

Fig. 1 Feynman diagrams contributing at tree- level to process (|T|). 
Fig. 2 Feynman diagrams contributing at tree- level to process (J|). 

Fig. 3 Cross sections of the processes ([!]) and (|3]) as a function of the Higgs mass, at y/s = 300 
GeV (upper plot), y/s = 500 GeV (central plot), y/s = 1000 GeV (lower plot), after the 
cuts p b rf > 1 GeV and cos 6^ < 0.95. 

Fig. 4 Cross sections of the processes (§) and (£|) as a function of the Higgs mass, at y/s = 300 
GeV (upper plot), y/s = 500 GeV (central plot), y/s = 1000 GeV (lower plot), after the 
cuts Pt' 7 > 1 GeV and cos 6^ < 0.95. 

Fig. 5 Distributions in invariant mass of the bb pair (upper left plot), in invariant mass of the 
667 system (upper right plot), in energy of the photon (lower left plot) and in missing 
transverse momentum (lower right plot) for processes © and (D, after the CUtS Prf > 1 
GeV and cos 9^ < 0.95, for a selection of Higgs masses: (a) at y/s = 300 GeV (bins 
of 2 GeV); (b) at ^ = 500 GeV (bins of 4 GeV); (c) at y/s = 1000 GeV (bins of 5 
GeV). In the first and third plot (clockwise) the upper histograms refer to rates from 
the non-radiative process (f|), whereas the lower histograms correspond to rates from the 
radiative process (p]). In the first plot the lowest order rates are shaded. In the fourth plot 
the upper histograms refer to rates from the complete set of diagrams describing process 
d5), whereas the lower histograms correspond to rates from contributions due to 'decay 
radiation' graphs only. 

Fig. 6 Distributions in invariant mass of the 4jet system (upper left plot), in invariant mass 
of the 4jet7 system (upper right plot), in energy of the photon (lower left plot) and in 
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missing transverse momentum (lower right plot) for processes (^) and (01), after the cuts 
Pt 1 > 1 GeV and cos6 l j 7 < 0.95, for various selections of Higgs masses: (a) at ^Js = 300 
GeV (bins of 2 GeV); (b) at s/s = 500 GeV (bins of 4 GeV); (c) at y/a = 1000 GeV (bins 
of 5 GeV). In the first and third plot (clockwise) the upper histograms refer to rates from 
the non-radiative process (^), whereas the lower histograms correspond to rates from the 
radiative process (§). In the first plot the lowest order rates are shaded. In the fourth plot 
the upper histograms refer to rates from the complete set of diagrams describing process 
(0), whereas the lower histograms correspond to rates from contributions due to 'decay 
radiation' graphs only. 
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Diagrams by MadGraph 
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e + e -> v v bb (y) 
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Fig. 5c 
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Fig. 6a 
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Fig. 6c 
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